Feeding compact objects: Accretion on all scales 
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Abstract. A major uncertainty in models for photoionised outflows in AGN is the distance of the 
gas to the central black hole. We present the results of a massive multiwavelength monitoring 
campaign on the bright Seyfert 1 galaxy Mrk 509 to constrain the location of the outflow 
components dominating the soft X-ray band. 

Mrk 509 was monitored by XMM-Newton, Integral, Chandra, HST/COS and Swift in 2009. 
We have studied the response of the photoionised gas to the changes in the ionising flux produced 
by the central regions. We were able to put tight constraints on the variability of the absorbers 
from day to year time scales. This allowed us to develop a model for the time-dependent pho- 
toionisation in this source. 

We find that the more highly ionised gas producing most X-ray line opacity is at least 5 pc 
away from the core; upper limits to the distance of various absorbing components range between 
20 pc up to a few kpc. The more lowly ionised gas producing most UV line opacity is at least 
100 pc away from the nucleus. 

These results point to an origin of the dominant, slow (u<1000 kms -1 ) outflow components 
in the NLR or torus-region of Mrk 509. We find that while the kinetic luminosity of the outflow 



2 



Jcllc Kaastra et al. 



is small, the mass carried away is likely larger than the 0.5 Solar mass per year accreting onto 
the black hole. 

We also determined the chemical composition of the outflow as well as valuable constraints 
on the different emission regions. We find for instance that the resolved component of the Fe-K 
line originates from a region 40-1000 gravitational radii from the black hole, and that the soft 
excess is produced by Comptonisation in a warm (0.2-f keV), optically thick (r ~10-20) corona 
near the inner part of the disk. 

Keywords. Active Galactic Nuclei, X-rays, outflows 



1. Introduction 

Active Galactic Nuclei (AGN) are powered by accretion onto a super-massive black 
hole. In addition to the accretion process, however, it has been found that many AGN 
contain strong outflows. The total mass loss associated with this outflow may even be 
larger than the mass that finally reaches the black hole, and the kinetic energy associated 
with the outflow may be large enough to regulate the growth of the black hole and tightly 
connected to it, the growth of the host galaxy. 

Different origins for these outflows or winds have been proposed, like accretion disk 
winds, thermal evaporation from the torus, or extended ionisation cones. These models 
all imply very different effects for the impact of the outflow on the environment, the 
feedback process. First, outflows launched close to the black hole (like accretion disk 
winds) have to overcome a larger gravitational potential to escape, hence should have 
much larger outflow velocities v, from thousands of kms -1 to significant fractions of the 
speed of light. Furthermore, outflows at larger distances r from the black hole (like e.g. 
torus winds) have a larger impact on the surroundings, because we have for the mass 
loss rate M: 

M/Q = N H m p rv, (1.1) 
with f2 the solid angle subtended by the outflow, m p the proton mass and Ah the total 
column density. The kinetic luminosity Lk is simply given by \Mv 2 . 

While column densities and outflow velocities can be obtained directly from the spec- 
trum, the distance should be obtained by other means. Given the ionisation parameter 
£ = L/nr 2 with L the ionising luminosity and n the hydrogen density, a known value 
for n immediately gives r, because L and £ are known from observations. Because the 
recombination time scale of the plasma scales with n _1 , the delay time of the ionisation 
state of the plasma with respect to variations of the continuum luminosity immediately 
yields the density and therefore the distance of the absorber. This is the principle that 
forms the basis for our large monitoring campaign on Mrk 509, a bright and luminous 
Seyfert 1 galaxy at z = 0.034. 

2. The monitoring campaign on Mrk 509 

We have observed Mrk 509 intensively during a 100-day campaign in 2009 (Kaastra 
et al. 2011a). The core of the campaign consisted of 10x60 ks observations with XMM- 
Newton, spaced 4 days, coincident with ten INTEGRAL observations of 120 ks each. This 
was followed by a 180 ks observation with the Chandra LETGS, simultaneous with 10 
orbits HST/COS. It was preceded with Swift monitoring (UV, X-ray) and supplemented 
with ground-based observations (WHT, Pairitel). Up to now, 12 refereed papers have 
been published on this campaign and more are under preparation. 

The lightcurve of Mrk 509 (Kaastra et al. 2011a) showed the expected variability with 
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a large (~60% increase) outburst in the middle of our campaign. Our first step was 
to determine the time-averaged spectrum. Because of the high statistical quality of the 
data, special analysis tools and refinements to the calibration needed to be developed, in 
particular for the RGS instruments of XMM- Newton (Kaastra et al. 2011b), but also for 
the EPIC camera (Ponti et al. 2012) and the Optical Monitor (Mehdipour et al. 2011) 
of XMM- Newton and the HST/COS spectrometer (Kriss et al. 2011). 

The time-averaged spectrum showed five different ionisation components producing 
significant absorption lines in the RGS spectrum (Detmers et al. 2011). These results 
were confirmed by the LETGS spectra (Ebrero et al. 2011). The quality of the data was 
high enough to prove for the first time that two of the ionisation components, C and D, are 
discrete structures with a FWHM of less than 35 and 80%, respectively. However, these 
components are not in pressure equilibrium. From component A to E, the gas pressure 
drops by two orders of magnitude. For a proper modeling of the ionisation structure 
the spectral energy distribution (SED) of the source needs to be known accurately, and 
thanks to our broad-band coverage from optical to hard X-ray wavelengths we were able 
to do this modeling. 

The next step was to constrain the distance r of the absorbers, by measuring the 
delayed response of the absorbers with respect to the continuum variations. To do this 
properly, we developed a time-dependent photo-ionisation model. Starting with equilib- 
rium, we followed the spectral variations over time and calculated the ion concentrations 
for different gas densities n. Comparing the predicted changes in absorber transmission 
for each of these models with the observed spectra, the distances could be constrained 
(Kaastra et al. 2012). For components C-E no significant variability during our campaign 
was found, leading to lower limits of their distance. However, component D showed ev- 
idence for variability on long time scales, by comparison with archival spectra taken in 
2000 and 2001, yielding a tight upper limit to its distance. For components C and E up- 
per limits were derived from the requirement that r < L/Nn£. Our measurements were 
not sensitive enough to constrain variability for components A and B, but Oiii is present 
in component A, and through direct imaging of its forbidden 5007 A line its distance 
has been determined to be ~3 kpc (Phillips et al. 1983). Furthermore, from variability 
constraints of the UV line velocity components, corresponding to ionisation components 
A and B, Arav et al. (2012) derived lower limits of ^200 pc for most of them. 

Fig. [1] summarises the distance limits. The dominant X-ray absorption components C- 
E originate from a region outside the inner torus radius. Thus, these outflow components 
are more likely explained by torus winds rather than accretion disk winds. 

Using several arguments, the total mass loss of the wind is between 0.26 and 2100 solar 
masses per year. This is large compared to the accretion rate of about 0.5 solar mass 
per year. However, due to the limited outflow velocities of the wind (up to several 100 
kms -1 , but not relativistic), the kinetic luminosity of the outflow is not very large. 

Our spectra also yielded the metal abundances of the outflow. Steenbrugge et al. (2011) 
measured the abundances of C, N, Ne, Mg, Si, S, Ca and Fe relative to iron with accuracies 
down to 8% for some elements. The ratio's agree fairly well with solar abundances, 
although we could not yet constrain the absolute abundances. We will determine those 
later using new HST/COS spectra covering better the Lyman series of hydrogen. 

The Galactic foreground absorption was studied by Pinto et al. (2012). The UV spec- 
trum shows seven discrete absorption structures. Some of these are associated to high- 
velocity clouds, others correspond to absorption within the spiral arms of our Galaxy. A 
detailed assessment of the role of dust and gas at various temperatures could be made. 

The UV lightcurve shows a distinct outburst near the centre of our monitoring cam- 
paign. The soft X-ray flux correlates very well with these variations in the UV flux, while 
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Figure 1. Distances of the absorption components in Mrk 509. Horizontal bars indicate the 
allowed ranges for each of the individual absorption components. The x-axis gives distances from 
the black hole, the y-axis the ionisation parameter £ for each component A-E. The constraint 
from the UV lines, a mixture of ionisation components A and B, is also shown. Vertically dashed 
lines indicate some radii of interest. 

the hard X-ray flux does not correlate at all with the other bands. This and the spectral 
shape led Mehdipour et al. (2011) to conclude that the soft X-ray excess in Mrk 509 is 
not caused by blurred reflection, but by comptonisation of the soft UV photons of the 
disk. Further detailed modeling by Petrucci et al. (2012) strenghtcned this picture. The 
scattering medium has a large optical depth of ^15, at a temperature of ~1 keV. 

Finally, Ponti et al. (2012) studied the Fe-K emission line, showing that the neutral 
component consists of a narrow, constant component, possibly originating from material 
between the outer broad-line region and the torus, and a broad variable component, 
originating from the inner broad line region. 



References 

Arav, N., Edmonds, D., Borguet, B., et al. 2012, A&A, 544, A33 
Detmers, R.G., Kaastra, J.S., Steenbrugge, K.C., et al. 2011, A&A, 534, A38 
Ebrero, J., Kriss, G.A., Kaastra, J.S., et al. 2011, A&A, 534, A40 
Kaastra, J.S., Petrucci, P.-O., Cappi, M., et al. 2011a, A&A, 534, A36 
Kaastra, J.S., de Vries, CP., Steenbrugge, K.C., et al. 2011b, A&A, 534, A37 
Kaastra, J.S., Detmers, R.G., Mehdipour, M., et al. 2012, A&A, 539, A117 
Kriss, G.A., Arav, N., Kaastra, J.S., et al. 2011, A&A, 534, A41 

Mehdipour, M., Branduardi-Raymont, G., Kaastra, J.S., et al. 2011, A&A, 534, A39 
Petrucci, P.-O., Paltani, S., Malzac, J., et al. 2012, A&A, in press (arXivl209.6438) 
Phillips, M.M., Bladwin, J. A., Atwood, B. & Carswell, R.F. 1983, ApJ, 274, 558 
Pinto, C, Kriss, G.A., Kaastra, J.S., et al. 2012, A&A, 541, A147 
Ponti, G., Cappi, M., Costantini, E., et al. 2012, A&A, in press (arXivl207.0831) 
Steenbrugge, K.C., Kaastra, J.S., Detmers, R.G., et al. 2011, A&A, 534, A42 



